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Abstract 

In supersymmetric theories, the main decay modes of scalar quarks are decays into 
quarks plus charginos or neutralinos, if the gluinos are heavy enough. We calculate the 
0{a s ) QCD corrections to these decay modes in the minimal supersymmetric extension 
of the Standard Model. In the case of scalar top and bottom quarks, where mixing 
effects can be important, these corrections can reach values of the order of a few ten 
percent. They can be either positive or negative and increase logarithmically with the 
gluino mass. For the scalar partners of light quarks, the corrections do not exceed in 
general the level of ten percent for gluino masses less than 1 TeV. 
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1. Introduction 



Supersymmetric theories (SUSY) [l], || are the best motivated extensions of the Standard 
Model (SM) of the electroweak and strong interactions. They provide an elegant way to 
stabilize the huge hierarchy between the Grand Unification or Planck scale and the Fermi 
scale, and its minimal version, the Minimal Supersymmetric Standard Model (MSSM) allows 
for a consistent unification of the gauge coupling constants and a natural solution of the Dark 
Matter problem [[|. 

Supersymmetry predicts the existence of a left- and right-handed scalar partner to each 
Standard Model (SM) quark. The current eigenstates, qi and qn, mix to give the mass 
eigenstates qi and q~2] the mixing angle is proportional to the quark mass and is therefore 
important only in the case of the third generation squarks In particular, due to the large 
value of the top mass m t , the mixing between the left- and right-handed scalar partners 
of the top quark, ti and t^, is very large and after diagonalization of the mass matrix, the 
lightest scalar top quark mass eigenstate t± can be much lighter than the top quark and all 
the scalar partners of the light quarks [[|. 

If the gluinos [the spin 1/2 superpartners of the gluons] are heavy enough, scalar quarks 
will mainly decay into quarks and charginos and/or neutralinos [mixtures of the SUSY 
partners of the electroweak gauge bosons and Higgs bosons] . These are in general tree-level 
two-body decays, except in the case of the lightest top squark which could decay into a 
charm quark and a neutralino through loop diagrams if the decay into a chargino and a 
bottom quark is not overwhelming ||. These decays have been extensively discussed in the 
Born approximation |J. In this paper we will extend these analyses by including the 0(a s ) 
corrections, which due to the relatively large value of the strong coupling constant, might 
be large and might affect significantly the decay rates and the branching ratiosQ 

The particular case of the QCD corrections to scalar quark decays into massless quarks 
and photinos has been discussed in Refs. P, [jj. In the general case that we will address here, 
there are three [related] features which complicate the analysis, the common denominator of 
all these features being the finite value of quark masses: (i) In the case of the decays of top 
and bottom squarks, one needs to take into account the finite value of the top quark mass 
in the phase space as well as in the loop diagrams, (ii) Scalar quark mixing will introduce 
a new parameter which will induce additional contributions; since the mixing angle appears 
in the Born approximation, it needs to be renormalized. (iii) The finite quark mass [which 
enters the coupling between scalar quarks, quarks and the neutralino/chargino states] needs 
also to be renormalized. 

The QCD corrections to the reaction q — > q\ analyzed in the present paper are very 
similar to the case of the reverse process, t — > tx° and t — > b\ + recently discussed in 
Ref. |TIJ (see also Ref. [|TT]]). During the preparation of this paper, we received a report by 



Kraml et al. |T2j, where a similar analysis has been conducted. Our analytical results agree 



1 If the gluinos are lighter than squarks, then squarks will mainly decay into quarks plus gluinos; the QCD 
corrections to these processes have been recently discussed in Refs.0, §]. 
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with those given in this paper[|. We extend their numerical analysis, which focused on the 
decay of the lightest top squark into the lightest charginos and neutralinos, by discussing 
the decays into the heavier charginos and neutralinos and by studying the case of bottom 
squarks and the SUSY partners of light squarks. 



2. Born Approximation 

In the Minimal Supersymmetric Standard Model |2|, there are two charginos xt = 1)2] 
and four neutralinos x° [i = 1-4]. Their masses and their couplings to squarks and quarks are 
given in terms of the Higgs-higgsino mass parameter /i, the ratio of the vacuum expectation 
values tan j3 of the two Higgs doublet MSSM fields needed to break the electroweak symmetry, 
and the wino mass parameter M 2 . The bino and gluino masses are related to the parameter 
M 2 [Mi ~ M 2 /2 and rrig ~ 3.5M 2 ] when the gaugino masses and the three coupling constants 
of SU(3)xSU(2)xU(l) are unified at the Grand Unification scale. 

The squark masses are given in terms of the parameters \x and tan /3, as well as the left 
and right-handed scalar masses Mq L and Mq R [which in general are taken to be equal] and 
the soft-SUSY breaking trilinear coupling A q . The top and bottom squark mass eigenstates, 
and their mixing angles, are determined by diagonalizing the following mass matrices 

KA 2_ (M? L +m 2 t +cos 2/3(1-14.) M 2 m t M t LR \ 

Mi ~{ m t M t LR Ml R + m 2 t + l cos 2f3 s 2 w Mf ) [L) 

U 2 _( M? L +ml + cos2/3(-§ + \s 2 w ) M 2 m b M b LR \ 

^i-y m b M b LR M? r + m 2 -\ cos 2(3 s 2 w M\ ) [Z) 

where M^ R in the off-diagonal terms read: M^ R = A t — fi cot (3 and M b LR = A b — fi tan f3. 

In the Born approximation, the partial widths for the decays tj — > fx®, t% — > hxj can be 
written as [q = t or b, and we drop the indices of the neutralino/chargino states] 



~ a 

ro(*i -> qx) 



4 mf 



(cl + 4 ) (m| - m\ - m 2 x ) - 4 c l L c l R m q m x e x 



X^imlmlm 2 ) (3) 



where X(x, y,z) = x 2 + y 2 + z 2 — 2 (xy + xz + yz) is the usual two-body phase space function 
and e x is the sign of the eigenvalue of the neutralino x- The couplings c l L R for the neutral 
current process, U — > tx°, are given by 



4 
4 






+ h | 




ce t J 




ce t 


\ + fn 




-s dt 


\ cet 



(4) 



2 We thank the Vienna group and in particular S. Kraml for their cooperation in resolving some discrep- 
ancies with some of the formulae and plots given in the early version of the paper Ref . (T^] . We also thank 
T. Plehn for checking independently the results. 
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and for the charged current process, ti — > &x H 



mi, £7. 



J2 



-Cflt 



svk M w cos /3 [ «6>i 



R 

-2 

R 



s w 



C0t 



mtV* 



J-2 



S0t 



V2s w M w sin/3 [ c 8t 



(6) 



In these equations, 6 t is the t mixing angle [which as discussed previously can be expressed 
in terms of the Higgs-higgsino SUSY mass parameter /i, tan (3 and the soft-SUSY breaking 
trilinear coupling A t ] with sg = sin9, eg = cos 8 etc.; s 2 ^ = 1 — c 2 ^ = sin 2 6 W and N,U/V 
are the diagonalizing matrices for the neutralino and chargino states [13| with 

N' jt = c w N n + s w N j2 , N' j2 = -s w N n + c w N j2 . (7) 

A similar expression eq. (3) can be obtained for the neutral and charged decays of bottom 
squarks, bi — > bXj and b — > txj 



r (&i -> qx) 



a 



4m~ i_ 

bi 



A^K, m 2 m 2 ) (8) 



with the couplings d L R in the neutral decay b — > bx° given by [9^ is the b mixing angle] 



R 

-2 
R 



brrib 

blTLf, 



s o b } 


+ h | 




c e b J 


~ s e b 




} + /« 


i s e b 


S01 


\ c b 



(9) 



b 

h 
Jr 



1 



N j3 



\pl M w cos (3 s w 
1 

'-1^1+*— M 



V2 
-V2 



CwSw 



N J2 



3 cw 



J-2 



and for the charged current process, bi — > tx 



ci 1 






c i J 


v^siyMpv sin/3 


{ } 


4 1 




m b U j2 


4 J 


% 1 J 


\/2 sw M w cos /? 



(10) 



(11) 
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In the case where the mass of the final quark and the squark mixing angle are neglected [as 
it is the case for the first and second generation squarks], the decay widths simplify to 



a 



Tofe -> qx) = ~r m m 1 



m 



2 \ 2 



ff 



'12) 



where the f^s [with now i = L,R since there is no squark mixing] in the case of the neutral 
decays, q — > qx°, are given in terms of the quark isospin l\ L and charge e g , by 



h 
Sr 



V2 
-y/2 



e q N'+(li L -e q s 



CwSw 



N> 2 



while for the charged decays, q — *> q'x + one has for up-type (down-type) squarks: 

h = Vji/sw {Uji/sw) , fn = 0- 



(13) 



(14) 



3. QCD corrections to Top Squark Decays 

The QCD corrections to the top squark decay width, eq. (3), consist of virtual corrections 
Figs.la-d, and real corrections with an additional gluon emitted off the initial t or final t 
[for the neutral decay] or b [for the charged decay] quark states, Fig. le. The 0(a s ) virtual 
contributions can be split into gluon and gluino exchange in the q—t~x [q = t, b] vertex as 
well as mixing diagrams and the t and t/b wave function renormalization constants. The 
renormalization of the q—t—x coupling is achieved by renormalizing the top/bottom quark 
masses and the t mixing angle. We will use the dimensional reduction scheme^ to regularize 
the ultraviolet divergencies, and a fictitious gluon mass A is introduced to regularize the 
infrared divergencies. 



3.1 Virtual Corrections 

The QCD virtual corrections to the U~x~q interaction vertex can be cast into the form 

6r=zef £ [<? jtL P L + (? jtR P R ] (15) 
i=g,g,mix,ct 

where G l , G~, G l mix and G l ct denote the gluon and gluino exchanges in the vertex, and the 
mixing and counterterm contributions, respectively. 

The contribution of the gluonic exchange [Fig. la] can be written as 

G l g>L>R = J LtR Fl + c\ L F l 2 (16) 



3 The quark mass and wave-function counterterms will be different in the dimensional regularization 
| |I4| and dimensional reduction schemes [ fi"5| . Since dimensional reduction is the scheme which preserves 
supersymmetry, we will present our results in this scheme. 
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with the form factors F[ 2 given by 



n 



B + 2m 2 q C -2 m| (C n - C 12 ) + 2 m\ C n 
-2 m, m x (C + Cn) 



(17) 



with q = t for the neutral and q = b for the charged decays; the two and three-point 
Passarino-Veltman functions, B = B (m 2 ,, A, mf.) and C. = C '..(m^, m|.,m^, m q , X 2 ,m~ ) 



can be found in Ref. 16 



The gluino exchange contributions [Fig. lb], are given by 

g I,l, r = -2£ <J(^f+44)*f t (44+44)^ 



fc=l,2 



- 44) F 6 fe =f (4*4 - v l4) F 7 k 



+4,L 



(«M + 44)*f =f (44 + v ¥i)F l i 



+{v\v\ - 44)Ff t (4*4 - v\c$Ff 



(1* 



with again g = £ for the neutral decay and g = 6 for the charged one; the form factors F{ k 7 
read 



iff 

tnfc 

^2,3 
rpik 
-^4,5 

cnfc 
^6,7 



m. 



m x [C + C i2 ] 



m x [±m, (C + Cn) + m t C 12 ] 
m 5 [m t C ± m q (Cn - C 12 )] 

ml C ± m t m q [C + C n - C 12 ] + [C n - C 12 ] + m y C 12 + B, 



(19) 



with the two- and three-point functions B = B (m 2 ,,m<j, m t ) and C. = C,,{m 2 q) m 2 ,m^, m~, 
m^,m 2 ). The couplings d k RL are given by 



jk 

a L,R 



c 



(20) 



for neutralinos, while for the charginos one has 

J «w 



4 

4 



J c e„ 

m t V j2 



m b U j2 



a/2 s w M w sin /3 \ sg b 



se 



a/2 s w M w cos (3 { ce, 



(21) 



The u~ and a~ couplings read 



4 = IH + ^J 



4 



\ K ~ s 9q ) 



(22) 



Finally, the mixing contributions due to the diagrams Fig. lc, yield the expressions 



G 



mix,L,R 



2 2 
m~ — ra- 
ti *2 



4ra t mj c 2 e t B (m { .,m t , m- g ) 
+ c 2 e t s 2 t (Ao(m|) - vVmfJ) 



(23) 



Therein, v4o is the Passarino-Veltman one-point function. Note that all these contributions 
are the same in both the dimensional reduction and dimensional regularization schemes. 

3.2 Counterterms 

The counterterm contributions in eq. (15) are due to the t and t/b wave function renormal- 
izations [Fig. Id] as well as the renormalization of the quark mass m t or m b and the mixing 
angle 8 t , which appear in the Born couplings. 

For the neutral decay process, U — > tx% the counterterm contribution is given by 



G 
G 



1,2 

ct,L 

1,2 
ct,R 



1 



c\ 2 (SZ R + SZ il2 ) + b{c 9t , -se t }5m t - bm t {s dt ,c et }59 t + f R {cg t , -s et }59 t 



Cr ( sz l + <5%. 2 ) + b W, c et }5m t + bm t {c 6t , -sg t }56 t - f L {s 9t , c 6t } . 



(24) 



whereas for the charged current process, ti — > bxt, one obtains, 



G 
G 



1,2 
ct,L 

1,2 
ct.R 



1 1, 

1 i, 



6Z b R + 6Z k .+2 



5m b 
m b 



+ 



rribUji 



- c 



R 



y/2 Sw Mw cos (3 
Suit V j2 



{se t ,c et }S9 t 



s w 



{s 9t ,Cg t }59 t 



\f2s w M w sin/3 
m t Vj2 



{se t ,c 6t } 



V2 sw M w sin (5 



{ce t ,-s et } S9 t 



(25) 



In the on-shell scheme, the quark and squark masses are defined as the poles of the propaga- 
tors and the wave-function renormalization constants follow from the residues at the poles; 



the corresponding counterterms are given by (see also Refs. [|T0| , |P2f 

5m n 



szi 



8Z~ t , 



£|(ra 2 J + ££(ra 



-^rH) 



mz 



SrK) + ^ / K) + 2E| / (mJ) 



£1K) + S «K) + 2S I'K) 



( s f)'K) 



(26) 



In the dimensional reduction scheme, the self-energies S and their derivatives up to a 
factor a s /3ii which has been factorized out, are given by 0, [L2|] 



mk 2 ) 



Vf)'(k 2 ) 



2B 1 (k 2 ,m q ,X) + (1 + c 2 e q )B 1 ( y k 2 ,mg,mq 1 ) + (1 - c 2Qq )B x {k 2 ,m g ,m q2 ) 



2B 1 (k 2 ,m q , A) + (1 - c 2dq )B 1 (k 2 ,m- g ,m qi ) + (1 + c 2 e q )B 1 (k 2 ,m g ,m 



<12J 



m 



4 B (k ,m q ,X)-\ s 2 e q (B (k , rrig, m qi ) — Bo(k } m~ g} m q2 )) 



m 



2 Bx(k\ m iv X) -2k B[(k 2 , m io A) + (mf +mj- k ) B' (k 2 , m t , m. 



B (k , m t , m-g) + 2 s 2e m t m g B' (k , m t , m s ) 



(27) 



Using dimensional regularization, the quark self-energies differ from the previous expressions 
by a constant; in terms of the their values in the dimensional reduction scheme, they are 
given by 



L,R Z i ^Sldim. reg. 



E% + 2 . 



(2f 



dim. reg 

Finally, we need a prescription to renormalize the t mixing angle 8 t . Following Ref. |L7 



we choose this condition in such a way that it cancels exactly the mixing contributions 
eq. (23) for the decay t 2 — > tx° 



56+ 



1 2 

m~ — mi 

H 12 



4m t mgC 2 e t B Q (m 2 2 ,m t ,mg) + c 2 g t s 2 e t (A (m 2 2 ) - A {m\)) . (29) 



Alternatively, since the lightest top squark t\ can be lighter than the top quark and 
then is more likely to be discovered first in the top decays t tiXo, one can choose the 
renormalization condition such that the mixing contributions are cancelled in the latter 
process; this leads to a counterterm similar to eq. (29) but with B (m~ , m t , m g ) replaced by 
B (m 2 ^,m t ,mg). The difference between the two renormalization conditions, 



Ase t 



4m t m-g c 2dt 

1 2 

m~ — m~ 

H *2 



BQ(m 2 v m t ,mg) - B (m 2 2 ,m t ,m~ g ] 



(30) 



is, however, very small numerically Indeed, if mi x is a few GeV away from m ( - 2 , one has 
9 t ~ — 7r / 4 and therefore c 2 g t ~ 0, leading to a difference which is less than one permille for the 
scenario of Figs. 2a/b. For degenerate top squarks, one has A59 = 4m t mgC 2 g t BQ(m 2 2 , m t , m§) 
which is also very small numerically [less than ~ 1% for the scenarios of Fig. 2.] 

The complete virtual corrections to the U — > qx decay width is then given by 



a a s 
6 m~ 7r 



Re \ (d L G\ + d R G R ) {mi -m\- m\) 



2 {c l L G R + c R G l L )m q m x e x 
8 



(31) 



The sum of all virtual contributions including the counterterms are ultraviolet finite as 
it should be, but they are still infrared divergent; the infrared divergencies will be cancelled 
after adding the real corrections. 



3.3 Real Corrections 

The contributions to the squark decay widths from the real corrections, with an additional 
gluon emitted from the initial t or final t/b states, can be cast into the form 



real 



2 a a s 
3m+. 7r 



C L ° l R m q m x 



e x [ ( m ii + m l~ m x ) J oi + m~. I 00 + m\ I n + I + I x ] 



+ (4 +c l R )[2(m^ + mJ-m|)(m|/ 00 +mJ/ii+Jo + /i) 



4 i^i _2 ^2^ T 



+2 [m q - (m x - m t 



I - /?] 



(32) 



where the phase space integrals I(rni,m q ,m x ) = / are given by |TJ 

„-2 
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re 1 
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4 m 2 m| 


hi = 


1 




4 m? 






+2Li 2 (l 
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1 




4 m? 
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4 m? 
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4 m? 
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4 m? 
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A x / 2 (m 
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A) 
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ln(/5 2 ) + 2 1n 2 (/3 2 )-ln 2 (/3 )-ln 2 (A) 
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.(33) 



m? — m 2 — m 2 + k m? — m 2 + m 2 — k m? + m 2 — m 2 — k 

Po — ^ j Pi — ^ ? Pi — — 



2 m q m x 



2 mi. m x 



2 mi m q 



(34) 



Our analytical results agree with the results obtained recently in Ref. [[LJ| . 
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4. QCD corrections to other squark decays 



4.1 Bottom Squark Decays 



In the case of the bottom squark decays, 6, — ► b\° and 6j — > £%~, the analytical expressions 
of the QCD corrections are just the same as in the previous section once the proper changes 
of the squark [mi. — > m bi ], the quark [q = b and q = t for the neutral and charged decays] 
masses and the mixing angles [6 t — ► 6 b ] are performed. The couplings for b decays are as 
given in section 2: for the d k L R couplings, one has in the case of the neutral decay b{ — > bx° 



a L,R 



(35) 



with c k L R of eq. (11), while in the charged decay 6j — > , they read 



4 

4 



4 
4 




J'2 



\/2s w M w sin/? 



C6 t 



V% s w M w cos P [ s dt 



(36) 



The counterterm contributions are the same as in eq. (24) with the change (t, t) — > (b, b) in 
the neutral decay; in the charged decay mode they are different due to different couplings 
(see also Refs. ITOl O): 
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ct.R 
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m t Vj 
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a/2 sw M w sin /? 
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\/2swMw cos/3 
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i s e b ,ce b } 



y/2swM w cos/3 
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(37) 



where again the c k L R are given by eq. (11). Except for very large values of tan /3, the b mixing 
angle [as well as the bottom quark mass] can be set to zero and the analytical expressions 
simplify considerably^. The case of the neutral decay b — > bx° is even simpler since one can 
also neglect the mass of the final b quark. In fact, the latter situation corresponds to the case 
of decays of first and second generation squarks into light quarks and charginos/neutralinos, 
which will be discussed now. 



4.2 Light Quark Partners Decays 

Neglecting the squark mixing angle as well as the mass of the final quarks, the virtual 
corrections of the processes q~i — * q\ [where the subscript % stands now for the chirality of 

4 In the absence of mixing, the left- and right-handed bottom squarks are, to a very good approximation, 
degenerate if Mq L — Mg R . In the rest of the discussion, Sl and bn [and a fortiori the partners of the light 
quarks and qn] will be considered as degenerate. 
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the squark, since in the absence of squark mixing one has q^R = qi ;2 \ are given by the sum 
of the gluon and gluino exchange vertices and the wave-function counterterm, plus the real 
correction. The total width can then be written as 



T = F 



4 ct 

1 + o- 1 (F g + F i + F ct + F r 

6 TV 



(38) 



where the decay width in the Born approximation T l has been given in eq. (12). In terms 
of the ratio k = m^/m^, the gluon exchange corrections are given by [A = 1/(4 — n) with n 
the space-time dimension, and \x is the renormalization scale] 



A 
2~ 



1 mi 



A 2 /m? 



A 2 /m 2 



= + 1 - In - ^ In 2 /J '"\ - In - Li 2 («) . 



(!-«)= 



1 - « 



The gluino exchange contribution, with 7 = m?/m~, is given by 



-ml -re 4 , 

K 1 — K 



1 [ 7 ln 7 -( 7 -l)ln( 7 -l)] + /t + 7 , 2 / 



l-« 



with 



(39) 



(40) 



m| (1 - k) 



Co(0,m?,m 2 ,m?,m?,0) . 



In terms of dilogarithms, the function 7 is given for kj < 1 by 



Lio 
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and for ivy > 1 one has 
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7 /t 
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(41) 



(42) 



7 



(43) 



The counterterm contribution, consisting of the sum of the squark and quark wave-function 
renormalization constants, reads 



ct 



A 
-— + 



7 



2 4(1 -7) 2 8 2 
1/ 2 ni 7-11 [ 2 7 -l 

"2 (7 - 1)111 ~r + i[(r^F 



/i 2 4 m~ 



+ 3 



In 7 . 



(44) 



Finally, the real corrections with massless quarks in the final state contribute 



-In 2 



A 2 /m 2 - 5, A 2 /m 2 - « (4 - 3k) 



q +-ln 

U-/t) 2 4 (1-k) 2 4(1 -k) 2 



In k 



-Li 2 («) - ln/dn(l - «) - 8 ^_^ - y +4 



(45) 
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We see explicitly that the ultraviolet divergences A/2 and the scale \i cancel when F % and 
F c \ are added, and that the infrared divergences ln 2 (A 2 /m~) and ln(A 2 /m~) disappear when 
F g , F ct and F r are summed. The gluino exchange contribution eq. (40) does not contain any 
ultraviolet or infrared divergences. The total correction in eq. (38) then reads 

Ftot = F g + Fg + ^ + F T 

1 /4 7 2 - 277 + 25 , 3k-5\ tt 2 nri/ ._, l /a lX 7 



+ H — 7, — 2 Li 2 (/c) — - (7 — 1) In 



8 \ 7 -l « - 1 / 3 v ' 2 w 7 

3 7 2 _ 4 +2 3 3k 2 -4k 1 

4(l- 7 )2 1117 ~ 2 ln(1 ~ K) + 4(«-l)2 " ln(1 ~ K) 

il n (l-K) + ^[7ln 7 -(7-l)ln(7-l)] + ^z^|/ 



+V K 7 



• (46) 



In the limit where the mass of the final neutralino or chargino is much smaller than the 
mass of the initial squark, the analytical expression of the QCD correction further simplifies: 

tot 4(7-l) 2 [ 2 W ; 7 4(7-1) 3 V ; 

Note the explicit logarithmic dependence on the gluino mass in the correction. This loga- 
rithmic behaviour, leading to a non-decoupling of the gluinos for very large masses, 

3 m 2 5 7r 2 

F tot = - In — § + - - — for m g > m~ q (48) 

4 m~ 2 3 

is due to the wave function renormalization and is a consequence of the breakdown of SUSY 
as discussed in Ref. ||. Had we chosen the MS scheme when renormalizing the squark/quark 
wave functions [i.e. subtracting only the poles and the related constants in the expression 
eq. (23)] we would have been left with contributions which increase linearly with the gluino 
mass. 

Our analytical results in the case of massless final quarks agree with the corresponding 
results obtained in Refs.[|, [/J, where the QCD corrections to the decay of a squark into a 
massless quark and a photino have been derived, after correcting the sign of F g in Ref. [Q; 
see also the discussion given in Ref. [0. 



5. Numerical Analysis and Discussion 

In the numerical analysis of the QCD corrections to squark decays, we will choose m t = 180 
GeV (consistent with ESI) and m& = 5 GeV for the top and bottom quark masses and a 
constant value for the strong coupling constant a s = 0.12 [the value of a s in the running 
from a scale of 0.1 to 1 TeV does not change significantly]; the other fixed input parameters 
are a = 1/128, M z = 91.187 GeV and s 2 w = 0.23 |2(J. For the SUSY parameters, we will 



take into account the experimental bounds from the Tevatron and LEP1.5 data [PHI , and in 
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some cases use the values favored by fits of the electroweak precision data from LEP1 W^. 

Fig. 2 shows the partial widths for the decays of the lightest top squark into the two 
charginos X1.2 an d a bottom quark [2a] and into the lightest neutralino Xi an d the sum of 
all neutralinos [the opening of the neutralino thresholds can be seen in the curves] and a top 
quark [2b]. In these figures, tan/3 is fixed to tan/5 = 1.6, a value favored by b-r Yukawa 
coupling unification [p3fl . The solid, dashed and dot-dashed curves correspond to the (M 2 , /i) 
values [in GeV]: (70, —500), (70, —70) and (300,-70) in Fig. 2a [which give approximately 
the same value for the lightest chargino mass, m + ~ 70 GeV] and (100, —500), (100, —100) 
and (250, —50) in Fig. 2b [giving an LSP mass of m x o ~ 50 GeV]. These values correspond 
to the scenarios M 2 <C M 2 — fi and M 2 3> and have been chosen to allow for a 
comparison with the numerical analysis given in [Q. The parameters in the t mass matrix 
are fixed by requiring m t - 2 = 600 GeV and varying Mr . The mixing angle is then completely 
fixed assuming M^ R = M~ th (0| « — 7r/4 except for mj very close to rn^ 2 ); in the bottom 
squark sector we have = 220 GeV, m~ b2 ~ 230 GeV and 6~ b ~ 0. 

Fig. 3 shows the magnitude of the QCD corrections relative to the Born width to the 
decays of the lightest top squark into charginos+bottom [3a/b] and neutralinos+stop [3c/d] 
for the scenarios described in Fig. 2a [for Figs.3a/b] and Fig. 2b [for Figs.3c/d]. For both 
the neutral and charged decays, the QCD corrections can be rather large and vary in a wide 
margin: from ~ ±10% for light top squarks up to ~ —40% for m% 1 ~ m ( - 2 and some (M 2 ,/i) 
values. 

The small spikes near mi x ~ 425 (530) GeV for x + 6 (x°t) decays are due to thresholds 
in the top squark wave function renormalization constants from the channel t\ gt. For 
the depicted range, this happens only for the value M 2 = 70 (100) GeV which leads to 
nig ~ 3.5M 2 ~ 245(350) GeV. Note, however, that when this occurs, the channel t\ — > gt 
becomes by far the main decay mode, and the chargino/neutralino modes are very rare. 

In Fig. 4 the variation of the QCD corrections for the decay ti — > bxf [4a] and ti — > tx\ 
[4b] is displayed as a function of the gluino mass, for two values of /i = —50 and —500 GeV 
and tan/3 = 1.6 and 20. The top squark masses are fixed to = 300 and = 600 GeV 
{9 1 = — 7r/4) and the b masses are as in Fig. 2. M 2 and hence the chargino and neutralino 
masses are fixed by m- g . The figure exhibits a slight dependence of the QCD correction on 
the gluino mass. For the chosen set of squark mass parameters, the variation of the QCD 
correction with /i is rather pronounced, while the variation with tan j3 is milder. 

Fig. 5 shows the partial decay widths for the decays of the lightest bottom squark [which 
in our convention is denoted by b\ and is almost left-handed] into the lightest chargino xi 
and a top quark [5a] and into the lightest neutralino Xi an d a bottom quark [5b]. As in 
Fig. 2, tan/? is fixed to tan/5 = 1.6 and = 600 GeV; the mass difference between the 
two squarks is ~ 10 GeV and we have for the mixing angle 9~ b ~ 0. The solid, dashed 
and dot-dashed curves correspond to the (M 2 , /i) values [in GeV]: (60, —500), (70, —60) and 
(300,-60) in Fig. 5a and (100, -500), (100, -100) and (250,-50) in Fig. 5b. The decay 
b\ — > txi is by far dominant when the channel b\ — > gb is closed, since its decay width is 
almost two orders of magnitude larger than the b\ — > LSP+bottom decay width. 
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Fig. 6 presents the magnitude of the relative QCD corrections to the decays b\ — > tx\ 
[6a] and b\ — > b\\ [6b] as a function of the bottom squark mass, for the same scenarios as 
in Fig. 5. Again, depending on the values of /j,, M 2 and , the QCD corrections vary from 
(±) a few percent up to —50%. 

Finally, Fig. 7 displays the QCD corrections to the decays of the SUSY partners of 
massless quarks into neutralinos, q — > qxo, as a function of the ratio k = m^/m"- for several 
values of the ratio 7 = m~/m~,7 = 1.2,1.5 and 2 [7a] and as a function of 7 for several 
values of k, k = 0.2, 0.5 and 0.8 [7b]. The quark mass and the squark mixing angle are set 
to zero and all squarks are taken to be degenerate. The corrections then depend only on the 
two parameters, k and 7 since the dependence on the other SUSY parameters factorizes in 
the Born term. The QCD corrections vary from small [most of the time negative] values for 
small k values and small gluino masses, up to ~ 20% near threshold. 

For the decays qi — > q'xf [the right-handed squark does not decay into charginos], the 
matrix elements in the chargino mass matrix do not factorize in the Born expressions and 
the QCD corrections further depend on the ratios Uji/Vji through the contribution Fg. This 
dependence is, however, rather mild since first the ratio Uj\jVj\ is of order unity in most of 
the relevant SUSY parameter space [in particular for > M 2 ] and second the contribution 
F g is small compared to the other contributions for gluino masses below 1 TeV. The QCD 
corrections for the decays q~L — > q'x are thus approximately the same as in the case of the 
decays into neutralinos. 

In conclusion: we have calculated the 0(a s ) QCD corrections to decay modes of scalar 
squarks into quarks plus charginos or neutralinos in the Minimal Supersymmetric Standard 
Model. We have paid a particular attention to the case of t [and also b] squarks, where 
mixing effects are important. In the case of top squark decays, the QCD corrections can 
reach values of the order of a few ten percent depending on the various SUSY parameters. 
They can be either positive or negative and increase logarithmically with the gluino mass. 
For the scalar partners of light quarks, the corrections do not exceed the level of ten to 
twenty percent for gluino masses less than 1 TeV. 
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Figure Captions 

Fig. 1: Feynman diagrams relevant for the 0(a s ) QCD corrections to the decay of a squark 
into a quark and a neutralino or chargino. 

Fig. 2: (a) Partial widths [in GeV] for the decays of the lightest top squark t\ into the two 
charginos xt an d xt an d a bottom quark; tan/5 = 1.6, m t ~ 2 = 600 GeV. The solid, 
dashed and dot-dashed curves correspond to the (M 2 ,/i) values [in GeV]: (70, —500), 
(70, —70) and (300, —70). (b) Partial widths [in GeV] for the decays of t\ into a top 
quark and the lightest neutralino Xi as well as the sum of all neutralinos [the thresholds 
can be read off the curves]; tan/3 = 1.6, mj 2 = 600 GeV. The solid, dashed and dot- 
dashed curves correspond to the (M 2 , lm) values [in GeV]: (100, —500), (100, —100) and 
(250,-50). 

Fig. 3: Relative size [in %] of the 0(a s ) QCD corrections to the decay rates (a) t± — > bxf, (b) 
ti — > bxt, (c) h — > txi and (d) ti — > J2i ^Xi, as a function of the top squark mass. The 
set of (M 2 , /i) parameters is as in Fig. 2a for Figs. 3a/b and as in Fig. 2b for Figs. 3c/d. 

Fig. 4: Relative size [in %] of the 0(a s ) QCD corrections to the decay rates (a) ti — > bxf 
and (b) ti — > tx\ as a function of the gluino mass. M 2 is fixed in terms of rrig with 
the GUT relation. The solid, dashed, dotted and dot-dashed curves correspond to the 
(tan/3,//) values [// is in GeV]: (1.6, -50), (20, -50), (1.6, -500) and (20, -500). 

Fig. 5: (a) Partial widths [in GeV] for the decays of bi into the lightest chargino xt an d a 
top quark; tan/3 = 1.6, m ( - 2 = 600 GeV. The solid, dashed and dot-dashed curves 
correspond to the (M 2 ,/i) values [in GeV]: (60, -500), (70, -60) and (300,-60). (b) 
Partial widths [in GeV] for the decays of b\ to lightest neutralino Xi an d a bottom 
quark; again tan/3 = 1.6, m t ~ 2 = 600 GeV. The solid, dashed and dot-dashed curves 
correspond to the (M 2 ,/x) values [in GeV]: (100, -500), (100, -100) and (250,-50). 

Fig. 6: Relative size [in %] of the 0(a s ) QCD corrections to the decay rates (a) bi — > txi 
and (b) bi — > &x? an d as a function of the bottom squark mass. The set of (M 2 , fx) 
parameters is as in Fig. 5a/b. 

Fig. 7: The size of the QCD corrections for the decays qi — > qx° as a function of the ra- 
tios m 2 x /m 2 q (a) and m?/m? (b). The solid, dashed and dot-dashed lines correspond 
respectively to 7 = 1.2, 1.5 and 2 in (a) and k = 0.2, 0.5 and 0.8 in (b). 
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(e) 



Fig. 1 
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